Feedback Control of Rabi Oscillations in Circuit QED 
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We consider the feedback stabilization of Rabi oscillations in a superconducting qubit which is 
coupled to a microwave readout cavity. The signal is readout by homodyne detection of the in- 
phase quadrature amplitude of the weak measurement output. By multiplying the time-delayed 
Rabi reference, one can extract the signal, with maximum signal-to-noise ratio, from the noise. We 
further track and stabilize the Rabi oscillations by using Lyapunov feedback control to properly 
adjust the input Rabi drives. Theoretical and simulation results illustrate the effectiveness of the 
proposed control law. 
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I. INTRODUCTION 

In control theory, the system to be controlled is com- 
pared to the desired reference, and the discrepancy is 
used to correct the control action^. In contrast to 
classical systems, where measurements do not alter the 
state of the system, quantum measurements will col- 
lapse the system instaneously into one of its eigenstates 
in a probabilistic manner: the "measurement-induced 
backaction" 2 . Although the quantum coherent feedback 
control has been proposed^ and extensively applied in 
quantum optics and cooling mechanical oscillators and 
so on^r— , the measurement-based feedback control still 
maintains a great interest. Based on the quantum tra- 
jectory theory, Wiseman and Milburn Ref. developed a 
quantum conditional stochastic master equation (SME) 
to describe the dynamics resulting from the feedback (of 
the measurement output at each instant) to the quantum 
system. SME has been a topic of considerable activity in 
recent years for it paves the way for studying real-time 
measurement-based feedback control 8 - - — in quantum in- 
formation processing and computation. 

Circuit quantum electrodynamics (i.e., circuit QED, 
where a superconducting qubit is coupled to a microwave- 
frequency resonator cavity; see, e.g., Ref. Il3 - [l5[ ) has 
been shown to be a promising quantum computing ar- 
chitecture. Circuit QED allows for rapid, repeated 
quantum nondemolition (QND) superconducting qubit 
measurement 2 *^ and also provides several simple high- 
fidelity readout mechanisms, such as using large measure- 
ment drive powers^, and using either quantum-limited 18 
or nonlinear bifurcation amplifiers^. Moreover, circuit 
QED is an excellent test-bed for implementing quantum 
feedback control in either the qubits or the microwave 
resonator— ~—. For example, a recent work 29 has been 
shown that quantum measurement-based feedback con- 
trol can reduce dephasing and remarkably prolong the 
Rabi oscillations. 

Here, we analytically derive a simple and 
experimentally-feasible measurement-based feedback 
control law for circuit QED to track and stabilize Rabi 



(a) 



nKOK^ 



I 



I 



Qubit 



^05h 



Readout 
Cavity 



Multiplier 




Time 

Delay 

T 






Rabi 
Reference 



(b) 

Rabi drive 

Read-out drive 



HWWWWWH 



Ons 



20 ns 



100 ns 



FIG. 1: (Color online) (a) Simplified circuit diagram of mea- 
surement and feedback control. A superconducting qubit (yel- 
low) is coupled to a microwave readout cavity (blue). The 
amplified output is homodyne-detected and the quadrature 
signal is then extracted from the noise by multiplying the 
time-delayed Rabi reference (green). The discrepancy is used 
to design the feedback control law to correct Rabi oscillations 
(red), (b) Schematic of the read-out drive to build up the 
photon population of the cavity and Rabi drive to stabilize 
the Rabi oscillation. 



oscillations. The paper is organized as follows. The next 
section contains a brief discussion of the circuit QED 
Hamiltonian, the quantum detection, and the stochastic 
master equation for the qubit. In Sec. Ill, we study the 
open-loop control of the Rabi oscillation. In Se. IV, we 
study the feedback control by the Lyapunov function 
method. We summarize our conclusions in Sec. V. 



II. CIRCUIT FOR MEASUREMENT AND 
FEEDBACK CONTROL 

As shown in Fig. 1(a), we consider a superconducting 
circuit QED system with a superconducting qubit cou- 
pled to a microwave readout cavity and driven by two 
external drives: (i) a read-out drive with amplitude e d (t) 
and frequency uj d near the cavity resonance frequency 
w c , and (ii) a Rabi drive with amplitude e r (t) and fre- 
quency uj r near the frequency of the qubit a;q t 19 ' 29 ~ — . 
The Hamiltonian of the entire system can be written as 



H = hui c a/a + h—<r z + hg(a/(T- + atj + ) 



+h[e d {t)e- Mdt a? + e^e^'V + h.c], (1) 

where a^ and a are the creation and annihilation op- 
erators for the microwave readout cavity, <r + and cr- 
ave the raising and lowering operators of the super- 
conducting qubit, and g is the coupling strength be- 
tween the cavity and the qubit. In the dispersive 
regime^, |A| = \co q — ui c \ >• g, by applying the dis- 
persive shift U — exp[g(acr + — aV_)/A], and moving 
to the rotating frames for both the qubit and cavity, 
U c = exp[— ia^auJdt], U q = exp[— ia z uj r t/2], with the 
rotating- wave approximation, the Hamiltonian in Eq. ([T]) 
becomes 



H c g = hA c a* a 

+h [e d {t)a) 



hxa [ aa z + n—o z + n—a x 



3(*)o] 



(2) 



where A c = cu c — ui d , X = 5 2 /A, Qr = 2e r {t)g / A and the 
Lamb-shifted qubit transition frequency uj q — Lj q —u) r +x- 
If the cavity state is coherent, and the microwave cavity 
decay rate is much larger than the qubit decay rate, k ^$> 
71 (that allows to decouple the qubit dynamics from the 
resonator adiabatically) , the state at time t is given by 
\g) <8> \a g (t)) or |e) <g> \a e (t)). Here \a g ^(t)) are coherent 
states of the cavity and, from Eq. ([2]), the field amplitudes 
are given by22, 



a g {t) 



-ie d (t) - i(A c - x)a g {t) - ~a g (t), 



a e (t) = -ie d (t) - i(A c + x)®e{t) - -^cx e (t). (3) 

Thus, these coherent states a g ( e ) act as "pointer states' — 
for the qubit. Based on homodyne detection, by applying 
the transformation 

P(t) = \e)(e\D[a e (t)} + \g)(g\D[a g (t)i 

with D[a] = exp[aa^ — a* a] as the displacement operator 
of the microwave cavity, the effective stochastic master 
equation for the qubit degrees of freedom is 

dp = -T~^% — [<r Zl p]dt- i— {<T x ,p]dt + 7i V [<j-\pdt 



70 + r d (i) 



Here 



and 



UJ ac (t) = U q + B(t), 



p{t) = o e (t) - a g (t) 



is the separation between the pointer states a g (t) and 
a e (t), r\ is the measurement efficiency, 7^ is the pure 
dephasing rate, T>[A] is the damping superoperator 



and 



Also, 



V[A]p = ApA^ - A^Ap/2 - pA^A/2, 



H[A]p = Ap+ ~ P A ] -{A + At) p. 



V[cr z ]pdt + ^Krj\/3(t)\H[cr z }pdWt. (4) 



T d (t) = 2 X Im[a s (*K(t)] 
is the measurement-induced dephasing and 

B(t) = 2 X Re[a g (t)a* e (t)} 

is the ac Stark shift. The innovation dWt is a Wiener 
process 7 with 

E [dW t ] = 0, and E[dW?} = dt. 

Due to the qubit decay 71 and dephasing 7^ + T d (t), the 
system must quickly lose its quantum features. 

A coherent drive is turned on for 20 ns to build up the 
photon population of the cavity and is then repeated ev- 
ery 100 ns (see Fig. 1(b)). The cavity pull is designed 
to be x/2tt = 5 MHz, and the cavity decay rate is 
k/2tt = 20 MHz. A homodyne detection of the readout 
cavity field, with the help of the distance /3(£) between 
the states \oi e (t)) and \a g (t)), can then be used to dis- 
tinguish the coherent states and thus readout the state 
of the qubit. By applying the P-transformation to the 
in-phase quadrature amplitude 

I = (ae-** + a)e 1 *) /2, 

with cf> the phase of the local oscillation, the homodyne 
measurement record coming from the microwave cavity 
becomes 

i(t) = y/Kfj\i3(t)\ (<j z (t))+m = s(t)+m, (5) 

where the qubit uncorrelated term y/Krj\p,(t)\sm[<fi + 
arctan(/i)], p = a g + a e , has been omitted. We have 
set the homodyne phase <j> to arg(/3), which corresponds 
to detecting the quadrature with the greatest separation 
of the pointer states. Here £(£) = dW t /dt is a Gaussian 
white noise, representing the shot noise, with spectral 
density Pf(u>) = 1. Usually, the quantum signal s(t) is 
very weak and the noise £(£) may be strong. The over- 
all objective is to make the system behave in a desired 



way by manipulating the input drive based on the mea- 
surement output. Here we expect to sustain the Rabi 
oscillations. To achieve this, the following steps are re- 
quired: Detect the signal s(t) from the noise £(£); re- 
construct X\{t) — Tr[cr x p(t)], X2(t) — Tr[a y p(t)], and 
xz{i) — Tr[a z p(t)], which are the three components of 
the Bloch vector for the ensemble qubit state based on 
the detected signal^; feedback the error signal between 
reconstructed state and the desired state, to design the 
feedback control law (the Rabi drive) thus minimizing 
the error. 



III. OPEN-LOOP CONTROL: NO FEEDBACK 

To see how the feedback Rabi drive will work, we first 
consider the open-loop control. Open-loop means that 
we do not use feedback to determine if the output has 
achieved the desired goal. One can simply drive the mi- 
crowave cavity with amplitude 

e r = 0, R S/2g 

to obtain the Rabi oscillation with Rabi frequency £Ir; 
but cannot correct any errors. To illustrate this, we 
have numerically simulated the microwave cavity field 
equation (3) and the superconducting qubit stochastic 
master equation (4) with the open-loop drive amplitude 
e r = fl R d/2g to obtain the expected Rabi frequency 
ttftj2Tt = 1 MHz, for four different measurement drives 
e d /2ir = 1 MHz, 4 MHz, 8 MHz, and 20 MHz. 

In Fig. 2, we show some of these numerical results for 
the open-loop control of Rabi oscillations with frequency 
VLr/2tt = 1 MHz. We set the initial state of the qubit as 
the excited state. Figure 2(a) shows the results averaged 
over 1000 realizations. In these results we set the mea- 
surement efficiency rj = 1, the qubit decay 7i/27r = 0.05 
MHz, and the pure dephasing rate r y<p/2ir — 0.1 MHz. 
The Rabi-drive amplitude e r = flud/2g and the fre- 
quency uj r = uj q + x, should be chosen carefully to make 
the Lamb-shifted qubit transition frequency equal zero. 
When acquiring information from the measurement, it 
of course induces significant backaction on the system. 
From Fig. 2(a), we see that for the small measurement- 
drive amplitude (ed/27r = 1 MHz, red solid curve), the 
qubit decays and pure dephasing dominates the evolu- 
tion. Thus, in this case, the measurement only causes 
small amplitude noise on the Rabi oscillation. How- 
ever, for the larger drive amplitude ed/2n — 4 MHz (not 
shown) and 8 MHz (black dotted curve) the measurement 
induces remarkable backaction on the qubit. 

We set td/2-K = 20 MHz to gain more insight into what 
is actually happening during the evolution of the Rabi 
oscillation with strong measurement-drive amplitude. As 
shown in Fig. 2, (a z ) exhibits decaying oscillations, in 
Fig. 2(a), when the drive is weak (ed/2ir = 1, and 8 
MHz) and discontinuous jumps between two levels, in 
Fig. 2(b), when the driving is strong (ed/2Tr = 20 MHz). 
Clearly, in the strong drive, the qubit will remain fixed 
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FIG. 2: (Color online) (a) ensemble-averaged behavior (over 
1,000 realizations) of the filtered signal (only plot (o" 2 (i)}) by 
continuous weak measurements of the open-loop controlled 
microwave readout cavity with Rabi frequency Qr/2tt = 1 
MHz. The Rabi-drive amplitude e r = Qn5/2g and the fre- 
quency uj r = ujq + x, which makes the Lamb-shifted qubit 
transition frequency equals zero. The qubit is initially in the 
excited state and the read-out drive amplitude are €d/2ir — 1 
MHz (red solid line), and 8 MHz (black dotted line), respec- 
tively, (b) shows the same quantity except just a single real- 
ization by strong measurement with measurement drive am- 
plitude of amplitude ed/2n = 20 MHz. 



in either z = +1 or —1. This is the Zeno effect. All 
these demonstrated that the open-loop control cannot 
compensate for the disturbances in the system. 



IV. FEEDBACK CONTROL 

We now propose a simple feedback control law allow- 
ing to compensate the dephasing of the superconduct- 
ing qubit, the measurement-induced backaction, and to 



maintain the coherence of the Rabi oscillations based on 
the above measurement scheme. The schematics of such 
feedback control is shown in Fig. 1. The amplified and 
filtered signal s(t) = (<r z (t)) is compared with the Rabi 
reference signal s*(t) — cos il R t, and the difference 



e(*) = s(i)-s*(t) 



(6) 



is used to generate the feedback signal u(i) that drives 
the microwave cavity in order to reduce the difference 
with the desired Rabi oscillations: e(t) — > (frequency 
tracking^). The difference e(t) evolves as 



i(t) = s(t)-s*(t) = E 



5<"<*» 



-8*(t) 



(7) 



= n R (t)(<T v (t))-'Yi(l + ((T z (t)))+n R smSl R t. 

Thus, we design the feedback control law (the Rabi-drive 
amplitude): 

«(*)= e(<) =- — (a y (t))- 1 [K x sign e(t) + K 2 e(t) 

--fi{l + (a z (t))) + n° R smn R t] , (8) 

where K%, K 2 > 0. Using the feedback-control law (jHJ in 
Eq. ©, we have 



i(t) = -Ki signe(i) - K 2 e{t). 



(9) 



Clearly, if e(t) > 0, then e(t) < 0; and if e(t) < 0, then 
e{t) > 0. 

The Lyapunov function metho d 36 ' 37 is usually em- 
ployed to prove the stability of an ordinary differential 
equation and widely used in stability and control theory. 
Here we can choose a simple Lyapunov function 



Obviously, 



u(t) = e\t)/2. 



v{t) > and u(t) = e(t)e(t) < 0. 



Then, the Lyapunov theorem tells us that every trajec- 
tory of Eq. © converges to zero: 



lim \s(t)-s*(t)\ ^0 



00, 



(10) 



which means the system is globally asymptotically stable. 
Now, the only problem is to choose K\ and K 2 . From the 
feedback control law in Eq. ©, we find that when s(t) is 
far from s*(t), a large K 2 is needed to make s(t) converge 
to s*(t) quickly. If s(t) is quite close to s*(t), signe(t) 
dominates the evolution, thus a small K\ is needed to 
reduce the error e(t). 

We have simulated the feedback loop designed above to 
maintain the Rabi oscillations with frequency £l° R /2ir = 
2.5 MHz. The measurement is set in the weak-driving 
regime, when the readout drive amplitude is e ( i/2Tr = 1 
MHz, where the measurement- induced backaction T^(t) 
and B(t) remain small. The control parameters K\ = 
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FIG. 3: (Color online) (a) Feedback-controlled ensemble- 
averaged (over 1,000 realizations) Rabi oscillations, which 
persist for much longer time than those with open-loop con- 
trol. The Rabi frequency Q.° r /2ty = 2.5 MHz and the read-out 
drive amplitude is ed/2n = 1 MHz. (b) Power spectral density 
for the averaged measurement of feedback-controlled Rabi os- 
cillations from (a) (red curve); the blue curve corresponds to 
the open-loop case with the same parameters of (a). 



5 x 10 6 and K 2 = 10 8 . The other parameters are the same 
as in the case of open-loop control. Figure 3(a) shows 
typical realizations of the feedback-controlled ensemble- 
averaged Rabi oscillations. Clearly, the feedback control 
can quickly track the reference Rabi signal and ideally 
fight against dephasing and the measurement-induced 
backaction. From Fig. 3 we can see that the feedback- 
controlled Rabi oscillations persist for much longer time 
than those with open-loop control. Finally, in Fig. 3(b), 
we compare the power spectral density of the averaged 
measurement record in feedback-controlled Rabi oscilla- 
tions (red curve) with the corresponding open-loop con- 
trol (blue curve). Both of them are centered at 2.5 
MHz. However, the feedback controlled spectrum has 
a needle-like peak at the Rabi reference frequency, while 
the open-loop controlled spectrum has a broad distribu- 
tion. Thus, we can precisely convert the amplitude of 
the Rabi microwave drive to a frequency. Clearly, the 
proposed feedback control has more advantages than the 
open-loop control, for stabilizing the Rabi oscillations in 
circuit QED. 



V. CONCLUSION 

In conclusion, we have proposed and analyzed a quan- 
tum feedback control method to stabilize the Rabi os- 
cillations in a superconducting qubit which is coupled 
to a microwave readout cavity. The control law can be 
conveniently tested in realistic quantum QED architec- 
tures. The output signal detection has been discussed 
and the maximum signal-to-noise ratio has been given. 



We have also analytically proven that the designed feed- 
back Rabi-drive amplitude can make the averaged fil- 
tered signal quickly converge to the reference Rabi signal. 
We have discussed the advantages of the quantum feed- 
back control, over the open-loop control, in stabilizing 
the Rabi oscillations. The proposed Lyapunov feedback 
control can be further applied to quantum state purifi- 
cation, quantum adaptive measurement, and quantum 
parameter estimation. 
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